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New physics in LHC Higgs boson pair production 

Matthew J. Dolan0 Christoph Englert0 and Michael Spannowskj0 

Institute for Particle Physics Phenomenology, Department of Physics, 
Durham University, Durham DHl 3LE, United Kingdom 

Multi-Higgs production provides a phenomenologically clear window to the electroweak symme- 
try breaking sector. We perform a comprehensive and comparative analysis of new electroweak 
physics effects in di-Higgs and di-Ifiggs+jet production. In particular, we discuss resonant di-Higgs 
phenomenology, which arises in the Higgs portal model and in the MSSM at small tan P, and non- 
resonant new physics contributions to di-Higgs production in models where the newly discovered 
Higgs candidate is interpreted as a pseudo-Nambu-Goldstone boson. We show that, for all these 
scenarios, a measurement of the di-Higgs and di-Higgs-|-jet final states provides an accessible and 
elaborate handle to understand electroweak symmetry breaking in great detail. 



I. INTRODUCTION 

Both ATLAS and CMS have observed a Standard 
Model-like Higgs boson [l[ at around 125 GeV Hi]. In 
the very same mass region, the combination of the D0 and 
CDF collaborations' data sets exhibits a SM-like Higgs 
excess with a local significance of 2.2(T The implica- 
tions of this newly-discovered particle have already been 
discussed in the context of the SM and beyond [Sl-Q • The 
combined local significance is mostly driven by an excess 
in the diphoton invariant mass, consistent with the SM 
Higgs boson within 2a. Therefore, we can expect that 
the observed particle bears some resemblance to the SM 
Higgs since — >■ /i — >■ 77 is sensitive to the special role 
of the Higgs particle in both the SM's gauge and Yukawa 
sectors and their interplay. Correlating this observation 
with electroweak precision data Q and measurements in 
the h -> ZZ, W~^W~ channels, which constrain the par- 
ticle's couplings to massive gauge bosons, we infer from 
fits to the data 0, 0] (most notably by the ATLAS them- 
selves 0) that the particle reproduces SM Higgs proper- 
ties within l-2tT. This agreement partially relies on bias- 
ing the fit towards the SM Higgs hypothesis byassuming 
a total decay width r{h — )• anything) ~ Vf^ Q and the 
absence of new degrees of freedom in gg, 77 — >■ h. These 
assumptions are, strictly speaking, neither theoretically 
nor experimentally motivated. A precise determination 
of the particle's couplings relaxing such assumptions is 
an LHC lifetime achievement, which will combine direct 
searches for heavy states that potentially run in produc- 
tion and decay loops and constraints of non-standard 
Higgs branching fractions. 

Deviations from the SM Higgs phenomenology even at 
the 10% level leave a lot of space for modifications of 
the Higgs sector by Beyond-the-SM (BSM) physics: new 
physics of roughly that size is largely unconstrained by 
the precise investigations of the SM at the Z mass pole. 
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Given that the corresponding BSM couplings need to be 
small, the current data does not provide constraints on 
weakly-coupled Higgs sector extensions beyond what we 
have already learned from LEP @. Currently, Monte 
Carlo-based analyses which target non-standard decays 
of the Higgs-like resonances [ifl-UJl suggest that branch- 
ing ratio limits of < 0(10%) can in principle be obtained 
at the LHC from direct measurements, depending on the 
characteristics of the non-standard decay. This bound 
might be too loose to efficiently probe interactions be- 
yond the SM. 

From this perspective, it is imperative to directly 
probe potential modifications of the electroweak sym- 
metry breaking sector, if phenomenologically possible, 
to fully exhaust the LHC's search potential to physics 
beyond the SM. One class of hadron collider processes 
which precisely serve this purpose is multi-Higgs produc- 
tion [Tsl . These processes are functions of the symmetry 
breaking potential's parameters and are, consequently, 
highly sensitive to the realization of electroweak sym- 
metry breaking. While triple Higgs production is be- 
yond the reach of the LHC experiments [l^, di-Higgs 
production can potentially be measured in rare decays 
pp — > /ift, — > 6677 [r^ . Only recently, the application 
of jet substructure techniques [ll] to di-Higgs produc- 
tion in boosted final states has uncovered sensitivity in 
pp — )■ hh{j) — > bbT'^T~ to both di-Higgs production 
and the trilinear Higgs coupling [l^. This approach is 
currently also being investigated by ATLAS [20| in the 
context of a LHC luminosity upgrade. 

Crucial to the findings of Ref. [l^ is accessing the 
small invariant di-Higgs mass phase space region which 
is mostly sensitive to the Higgs trilinear coupling with 
moderately boosted Higgses px ~ 100 GeV. The sensi- 
tivity can be augmented by accessing collinear di-Higgs 
configurations by recoiling the di-Higgs system against a 
hard jet This configuration is extremely sensitive 

to modifications of the trilinear Higgs coupling since it 
does not suffer from the kinematical shortcomings that 
are present in the inclusive di-Higgs final state, where the 
Higgs particles are produced back-to-back. Promising re- 
sults to measure the di-Higgs cross section have also been 
found for extremely boosted bbW^W~ production [2]| . 
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Motivated by the recently-unravelled sensitivity to di- 
Higgs production at the LHC, we perform a comprehen- 
sive and comparative analysis of new physics interactions 
in LHC di-Higgs and di-Higgs-|-jet production in this pa- 
per. We divide our discussion into two parts. We discuss 
resonant di-Higgs(-t-jet) signatures in Sec. Ill A| where we 
first analyze a simple extension of the Higgs sector via 
the so-called Higgs portal [22[. We subsequently discuss 
prospects to constrain the MSSM Higgs sector at low 
tan /? via resonant production of a heavy Higgs H decay- 
ing to hh. 

In Sec. IHII we discuss the phenomenology of non- 
resonant new physics contributions to di-Higgs produc- 
tion in composite Higgs and dilaton models (to make this 
work self-contained we briefly introduce the basics before 
we comment on the phenomenology). This broad class 
of pscudo-Nambu-Goldstone theories comprises many in- 
teresting features in a phenomenologically well-defined 
framework. Both these models introduce new degrees 
of freedom to di-Higgs and di-Higgs-|-jet production and 
modified trilinear couplings compared to the SM, while 
the composite Higgs scenarios also introduce new tihh 
interactions. Comparing these models to the SM ex- 
pectation provides a consistent framework to constrain 
the electroweak symmetry breaking potential with future 
measurements at the y/s = 14 TcV LHC. 

Throughout this paper, we produce events and lead- 
ing order cross sections using an in-house Monte Carlo 
code that is based on the Vbfnlo HI] and Fey- 
nArts/FormCalc/LoopTools [11] frameworks. 



II. RESONANT NEW PHYSICS: FROM THE 
HIGGS PORTAL TO SUPERSYMMETRY 

A. Di-Higgs production in the Higgs portal 
scenario 

The Higgs portal scenario [g^l is a convenient and theo- 
retically consistent way to generalize the SM in its mostly 
unconstrained parameters (such as Higgs boson's total 
and hidden decay width) in a minimal approach [25| . Re- 
alizing that 'I'^'I's transforms as a gauge singlet, where 
$5 is the SM Higgs doublet, there is a plethora of SM ex- 
tensions with hi ghly modified and interesting LHC phe- 



nomenology [10|, lUj, |26| . In a 'mirrored' approach [27 



the Higgs portal potential reads 

V = mII^sP + As|$s|' + 1^h\^h? + Ah|$h|^ 

+ t^^\<^s?\^h\\ (2.1) 

where we have introduced a hidden sector Higgs field 
<^H- The Higgs portal model allows to identify a viable 
dark matter candidate in the hidden sector [28|, whose 
potential LHC phenomenology has been explored in [2^ . 

After symmetry breaking, which is triggered by the 
Higgs fields acquiring vacuum expectation values (vevs) 
I^'s.hI = vs^h/V2, the would-be-Nambu-Goldstone 
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FIG. 1: Mass of the heavy Higgs state H il mu ^ 125 GeV, 
and consistency with 5, T, U [30|] , unitarity and current AT- 
LAS/CMS results is imposed. The density of the model points 
must not be interpreted as a probability measure. 



bosons are eaten by the , Z fields and the correspond- 
ing directions in the hidden gauge sector, and the only 
effect (in unitary gauge) is a two-dimensional isomctry 
which mixes the visible and the hidden Higgs bosons: 



h = cos xHs + sin x Hh 
H = — sin X Hg + cos x Hh , 



(2.2) 



where x is a function of the portal potential parameters 
Eq. (|2.ip (for details see, e.g., [1^). For the remainder 
of this section we choose niH > ruh = 125 GeV. 

Electroweak precision measurements and unitarity re- 
quirements of longitudinal gauge boson scattering and 
massive quark annihilation to longitudinal gauge bosons 
suggest that, if such a model is realized in nature, 
then the mixing should preferably be far from maximal, 
cosx^ ~ 1, which for generic perturbative choices of the 
potential XsjXviVx ^ '^''^ results in a typically small 
mass splitting between the physical Higgs states h, H. 
Admitting some tuning, a larger mass splitting can be ar- 
ranged, which results in a clean LHC phenomenology of 
narrow trans-TcV resonances [sij . Small mass splittings 
imply a phenomenologically more involved situation since 
the light Higgs bosons are produced with small transverse 
momentum in di-Higgs production. Nonetheless, given 
the vastly enriched Higgs sector phenomenology, we can 
still study the Higgs portal in sufficient detail to fully 
reconstruct the Higgs potential Eq. (|2.ip [25j . Crucial 
in this reconstruction algorithm is the measurement of 
the invisible Higgs decay branching ratio [l^, HI] . It can 
be immensely improved by a possible observation of a 
cascade decay H — > hh. Additional information from ob- 
serving all multi-Higgs signatures (and the trilinear cou- 
plings especially), if phenomenologically accessible, can 
be used to further constrain or even rule out the simple 
portal extension. 

Expanding Eq. (j2.ip around the vacuum expectation 
values, we get the trilinear couplings relevant for di-Higgs 
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(a) (T/iT(portal) and invariant di-Higgs 
mass distribution for pp — > hh + X at the 
LHC 14 TcV. 
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(b) (T/(T(portal) and invariant di-Higgs 
mass distribution for pp hH -|- X at the 
LHC 14 TeV. 
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(c) (T/(T(portal) and invariant di-Higgs 
mass distribution for pp — > HH -|- X at 
the LHC 14 TcV. 



FIG. 2: Upper panels: cross sections in the portal scenario for the parameter point mentioned in the text. We scan over the 
multiples of the trilinear couplings Eq. (|2.3p that are in one-to-one correspondence with diagrams involving the h, H propagators 
and show contours relative to the central expectation Eq. H2.4[) . Lower panels: invariant di-Higgs mass distributions. 



production* : 

hhh : ■i/2{2\Hs\vH + 2Xsc\vs 

+ Vx^x^x'^H + ij^c^slvs) , (2.3a) 
HHH : 3/2(2Ah4i,h - 2Xsslvs 

+ VxCxslvH - r]^cls^vs) , (2.3b) 
hHH : 2{3Xh - riy^)c\s^Vh + 2(3A5 - rj^y^s^^vs 

+ Vx^x""" + Vx^x'"s > (2-3c) 
hhH : 2{3Xh - riy^)c^s^^VH - 2(3A5 - ri^)c\s^vs 

+ -nx'^\vH --qx^^vs , (2.3d) 

where c^^ = cosx and s^^ = sinx- Current observations 
leave open a lot of parameter space for such signatures 
to be relevant at the LHC. In Fig. [T] we scan over the pa- 
rameters of the Higgs portal potential enforcing unitarity 
and electroweak precision constraints, as well as current 
limits from the ATLAS and CMS experiments [IH. If 



the heavier Higgs mass is mn > 250 GeV, there are pa- 
rameter choices such that the sin^ x suppression of the H 
decay to SM matter from Eq. (|2.2p renders the prompt 
decay of H to observable SM matter subdominant to the 
cascade decay H — > hh. This can be traced back to large 
trilinear couplings 0{vh,vs) that arise as a consequence 
of electroweak symmetry breaking. Therefore, there is 
the possibility to constrain the portal model by measur- 
ing the trilinear couplings in resonant and non-resonant 
pp — > hh, hH, HH + X production. 

In Fig. [21 we show a scan over the cross sections 
of pp — ?> hh, hH, HH -> (visible) as functions of the 
involved trilinear couplings for a exemplary parameter 
point vs ^ 246 GeV, vh ^ 24 GeV, ruh = 125 GeV, and 
TTiH ■^^^ 255 GeV, Th = 24 GeV. The central inclusive 
cross section values at leading order implying (prompt) 
visible final states are 

pp hh + X : 44.4 fb (2.4a) 
pp^Hh + X : 5.57 fb (2.4b) 
pp^ HH + X : 667 ab (2.4c) 



*Triple Higgs production, which is sensitive to the modified Higgs 
quartic couplings yields phenomenologically irrelevant cross sec- 
tions just like in the SM [3. 



(the SM cross section is 16 fb). Comparing to the 
NLO QCD corrections in the context of the (MS)SM 
by running HiGLU [1^ and HPAIR 34 [, we can expect 
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FIG. 3: Invariant mass distribution for pp ^ hh -\- j -\- X m 
the portal scenario. 



an enhancement of the cross section by about K = 

^NL0/^L0^2. 

For pp hh + j + X with pTj > 80 GeV we calcu- 
late a leading-order cross section of cr = 10.1 fb (Fig. [3]) 
which should be contrasted to a SM leading-order cross 
section of cr = 2.58 fb, which can be isolated from the 
background [I^. Hence, the measurement of the one 
jet-inclusive cross section will assist in formulating con- 
straints on such a model. 

Note that, pp HH + X ^ visible is naively sup- 
pressed ^ sin^ X- Therefore, for the bulk of the por- 
tal parameter space, heavy di-Higgs production (and di- 
Higgs-|-jet production different from pp — >■ hh + j + X) is 
phenomenologically inaccessible at too small rates, with 
no space left for kinematical signal-over-background S/ B 
improvements. 

Summary: The Higgs portal scenario offers the pos- 
sibility of large enhancements in the di-Higgs produc- 
tion rate, from both resonant and non-resonant (via 
changes in Xhhh) new physics. Extracting the rate for 
pp ^ h* ^ hh using the boosted kinematical techniques 
from our previous paper [Toj along with measuring the 
resonant peak in the di-Higgs invariant mass spectrum 
will aid in the full reconstruction of the Higgs portal la- 
grangian by correlating these two independent measure- 
ments. This strategy is assisted by the cross section's 
large dependence on Xhhh ■ A high luminosity analysis of 
hh and hh + j production can also facilitate a measure- 
ment of the trilinear coupling in this model. 



B. The MSSM at small tan /3 



The Higgs portal model of Sec. IH Al bears some re- 
semblance to a generic two Higgs doublet model, and 
therefore our findings are also relevant for searches for 
supersymmctry in the context of the MSSM and its ex- 
tensions. 

The trilinear couplings of the Higgs bosons in the 



Xhhh = 
Xuhh = 
XuHh — 



3 cos 2a sin(/? + a) 

2 sin 2q! sin(/3 -I- a) — cos 2a cos(/3 -I- a) 
- 2 sin 2a cos(/3 -|- a) — cos 2a sin(/3 -|- a). 



(2.5) 



up to radiative corrections, details of which can be found 
in the second reference of [l5[, tan (3 = Vu/ Vd the ratio 
of vevs of the two MSSM Higgs doublets, and a diag- 
onalizes the Higgs mixing matrix. The above couplings 
are in units of Aq = M^jv. In principle, disentangling 
the contributions proportional to Xuhh and Xhhh in dou- 
ble Higgs production would allow a reconstruction of the 
angles a and /3. 

We observe that when /? is small and we are near the 
decoupling limit where a ~ /3 — 7r/2 then the Xuhh is 
proportional to cos/3. Thus when Ivrih < mu < 2mt 
H has a large branching ratio into a pair of Higgses hh, 
similar to the Higgs portal model in Sec. Ill Al Probing 
the dihiggs final states is thus probably the best way 
of finding H if tan /3 is low. The presence of squarks 
can further enhanced the production by running in the 
loops^. 

Achieving a Higgs mass of 125 GeV at such low values 
of tan /3 requires exceptionally heavy stop masses and 
mixings. Scanning over the squark masses and mixings, 
we find that > 50 TcV in order to achieve ruh ^ 
125 GeV. These spectra are characteristic of 'mini-split' 
SUSY, which has recently been advocated in [s^ , which 
suggests that the weak scale is tuned and supersymmctry 
present at higher energies. However, it is unusual to have 
all the scalars heavy except the extra Higgses. This would 
require the presence of a cancellation between rnj^^ and 
mfj^ if these quantities are large like the other scalar soft 
terms, or else that they have some suppression relative 
to the squark and slepton masses. 

Moving beyond the MSSM, another possibility is that 
tan /? is low and that a large contribution to the mass of 
the lightest (SM-like) Higgs boson comes from an extra 
singlet field S with superpotential couplings XSH^Hd. 
This induces an extra contribution to the Higgs mass 
cx sin^ 2/3 which is enhanced at low tan /3; this is the 



so-called A-SUSY scenario of the NMSSM [Sg. We fo- 
cus exclusively on the MSSM here, however wc expect 
the phenomenology to be similar in the NMSSM if the 
singlet-like states are heavier than the MSSM Higgses. 

We find a point with tan/3 = 3, and adjust the scalar 
masses until we achieve mh ^ 125 GcV. Wc set the 
mass of the other CP-even boson of the MSSM H to be 
290 GeV. In this regime the branching fraction BR(i/ ^• 
hh) ~ 45%, and the decay width is Th = 0.25 GcV. 
The other main partial decays are into hh (12%), W^W~ 



^Note that, depending on the color charge assignment di-Higgs pro- 
duction can be enhanced compared to single Higgs production |35|l . 
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FIG. 4: Invariant mass distribution of the (a) hh and the (b) 
Hh system for MSSM-hke production at low tanp. For details 
see text. 



(28%) and ZZ (12%). We could increase the branching 
ratio into two Higgses further by decreasing tan (3, at the 
cost of increasing the scalar masses. Using a suitably 
modified version of Vbfnlo we find the leading order 
production cross-section a{pp H ^ hh) = 246 fb. We 
also calculate the cross-section for (j{pp — ^ iJ — >■ Hh). 
This is suppressed by the off-shell H in the s-channel, 
and by the fact that the XuHh coupling is suppressed 
relative to the Xnhh coupling. We find the cross-section 
for this process to be 4.5 fb, too low for observation given 
h has SM-Higgs-like branching ratios. 

We can separate the large contribution H — > hh by 
reconstructing the di-Higgs invariant mass which exhibits 
a peak at mn- This allows us to extract the cross-section 
for pp H ^ hh, and after cutting around the peak the 
remainder of the events are due to pp ^ h ^ hh. As 
in the Higgs portal model, this process can be extracted 
using the techniques from our previous paper, allowing 
constraints to be put on a and /3. The invariant mass 
distribution and rate for the hh + j final state are also 
similar to the portal scenario. Fig. [3] 

Summary: The di-Higgs phenomenology in the IVISSIVI 
at low tan /3 is similar in many respects to that of the 



Higgs portal model. JVIeasurements of the resonant and 
non-resonant contributions to di-Higgs production allows 
a reconstruction of the parameters a and /?. 



III. NONRESONANT NEW PHYSICS: 
PSEUDO-NAMBU-GOLDSTONEISM 

Apart from softly-broken supersymmetry, strong in- 
teractions are the only other constructions which can 
cure the naturalness problem (if only partially) with phe- 
nomenologically testable implications. 

A well-known example of electrowcak symmetry break- 
ing from strong interactions is technicolor (TC) where 
iTiw ~ / where / is the "pion" decay constant. The 
techni-E and techni-p resonances will have masses of the 
order of the TC confining scale, which can be much larger 
than the electroweak scale, Atc ^ /■ This usually trig- 
gers a tension with curing the quadratic energy diver- 
gence in perturbativc longitudinal gauge boson scatter- 
ing, which demands at least a single light degree of free- 
dom. An illustrative example which incorporates such 
a state is easily constructed from the holographic inter- 
pretation of a bulk gauge theory broken by boundary 
conditions in a Randall-Sundrum background [ssj^: The 
appearance of the infrared brane signals the spontaneous 
breakdown of conformal invariance in the dual picture 
poj . This is accompanied by higgsing of a symmetry, 
which is weakly gauged into the strongly-interacting sec- 
tor. On the one hand, such a "higgsless" theory does not 
have light scalar degrees of freedom analogous to the SIVI 
Higgs boson. On the other hand, stabilizing the compact- 
ification moduli via the Goldberger-Wise mechanism [4l| 
lifts the zero mass radion, which couples to the conformal 
anomaly 



^^z,z^ 

cos"^ 



f 



(3.1) 



In the CFT picture wc identify a pseudo-dilaton, which 
has an impressive resemblance to the SIVI Higgs boson as 
a consequence of its couplings. In this sense, the dilaton 
mimics a light Higgs boson because the mass terms are 
the source of scaling violation. 

Difl:crent to this approach is the interpretation of the 
entire Higgs multiplet as a set of Nambu-Goldstone 
bosons. There are multiple ways to construct such a 
model consistently, ranging from collective syrnmetry 
breaking 42 1 to holographic Higgs models [i^, |43| which 



vary in their details and symmetry content. Common to 



■t Owing to the large N and large 't Hooft coupling limit [3ij| of 
AdS/CFT, it is intrinsically difficult to construct a fully reaUstic 
model in terms of electroweak precision measurements. 
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all these realizations is the breaking of a global symmetry 
pattern by gauging a subgroup of the strongly interacting 
sector. 

While there are parameter choices for both scenarios 
which are consistent with the SM in their single Higgs 
phenomenology, the measurement of the di-Higgs(+jet) 
production can be a key discriminator between these dif- 
ferent non-resonant realizations. 



A. Di-dilaton production 

We first discuss the implications of interpreting the 125 
GeV boson as a pseudo-dilaton [4^|4^. We note that 
there is a substantial number of options in modelling the 
electroweak sector using strong interactions, and thus the 
conclusions of this section should be taken as illustrative 
rather than definitive for this class of models. 

The pseudo-dilaton is associated with the spontaneous 
breaking of scale symmetry at an unknown scale /, and 
we denote this field by x- The couplings of the pseudo- 
dilaton to massive Standard Model particles are deter- 
mined by its coupling to the trace of the SM energy 
momentum-tensor T^i,, Eq. p.ip . The couplings of the 
dilaton to the massive SM particles are thus the same 
as those of the SM Higgs, but rescaled by a factor of 
v/ f. The couplings of the pseudo-dilaton to gluons and 
photons are given by 



C 



X,masslcss 



CtEAI 
87Tf 



as 



(3.2) 

where cem.g arc anomaly coefficients. The precise value 
these take depends on what further assumptions are 
made about the UV dynamics of the theory and what 
heavy colored and electromagnetically charged states are 
present. We assume that the dilaton couples to photons 
and gluons via the fuh QCD/EM beta- function |47ll48l|. 
We also consider the model of [i^ which studies the same 
system with an extra family of quarks. In both these 
cases we can find SM-like behavior with an enhanced 
a X BR into photon pairs. 

In the fully composite scenario, the dilaton couples to 
massless gauge bosons via the full beta-function, we have 
cem = —17/9 and cs = 11 — 2n//3 where nj = 5 [47| is 
the number of light quarks. In the four-family model 
we have cem = —6/5 and Cs — 4/3, and obtain a 
similar single-dilaton phenomenology. The production 
cross-section of the dilaton can be enhanced by orders 
of magnitude relative to the SM value. However, as 
the dominant decay channel then becomes X ~^ 99^ the 
cross-section times branching ratio of the observable fi- 
nal states X — ^ //j X — ^ W and X ~^ 77 can still 
be close to their SM values, depending on the scale / 
(following arguments similar to the ones presented in 
Ref. [iSlsSEl)- 

There will also generally be dimension six operators. 
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FIG. 5: Dilaton production from gluon fusion with current 
limits of the /177 coupling analysis [g] included. 



the most interesting of which is 



C 



a. 



4^/ 



2CxxGGX^ {Gl^yf 



(3.3) 



We define the DQ operator with a minus sign, so 
that c^xGG > complies with the low-energy effective 
Higgs theorems [13, [13, 1131 paradigm: Integrating out 
the heavy top quark, we obtain an effective interaction 
C - C^^G"'''' log(l + h/v) in the SM. 

It is important to keep in mind that the higher di- 
mensional interactions with the gluon and the photon 
fields arise from integrating out the conformal dynamics 
and need not follow the LET paradigm, which predicts a 
unique coupling structure of the //."G^^ interactions as a 
consequence of m cx {h) for all fundamental masses in the 
SM§. We also explore the possibility that the dimension 
six operator is negligible, by setting c^xGG ~ 0- 

For the fully composite model we find that for / = 
850 GeV that a x BR of the dilaton into massive final 
states are very similar to those in the Standard Model, 
and X — >■ 77 it is 1.55 times the Standard Model value. 
For gluons we find a x BR is enhanced by a factor of 
approximately 150. This agrees with values obtained 
from recent fits of experimental data in |45| . We 
show in Fig. [5] the a x BR for massive states and for 
77, normalized to the SM values . We also include a 
blue horizontal band indicating the signal strength in the 
diphoton channels from combining the ATLAS and CMS 
searches 0, 3 • In the four- family case we obtain similar 
resuhs for / - 570 GeV. 

The dilaton's total decay width is approximately 
5 MeV and very similar to the SM. Upper limits on the 
Higgs width are difficult to assess experimentally [56[ and 
will eventually be limited by large systematic uncertain- 
ties [57| . Constraints on the dilaton model arising from 



^It is intriguing to realize that there is, in fact, a connection between 
LET and the vanishing trace anomaly Eq. l|3.1|l for infrared photons 
limQ2^0 (011^/^177) = mil. 
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such measurements will be too loose to rule this model 
out. As we will see, investigating multi-dilaton produc- 
tion provides the missing handle to constrain the model 
when consistency with single Higgs observations prevails. 

We introduce explicit sources of scale symmetry break- 
ing [131 through the operator XoO{x), where the scal- 
ing dimension of the operator O ^ A induces non- 
derivative trilinear interactions for x- When the operator 
O is nearly marginal, so that its anomalous dimension 
7 = \Aq — 4| ^ 1 and one writes the trilinear coupling 

as ■|-^x'^j one obtains for A 

A = (Ao + 1) + . . . , Ao < 1 (3.4) 

where we must have at A > 2 by the conformal alge- 
bra and unitarity. If Ao = 2 we obtain the Standard 
Model result, rescaled by the ubiquitous factor of v/ f. 
Another possibility is when 7^1 when one obtains 
A = 5, 66% larger than the SM trilinear up to factors of 
v/ f. There are also interesting anomalous four-derivative 
interactions in the low-energy dilaton theory [s^ [s^ , 

Z:^^^™ D 2{auv - aiR){2{dx)^ax - (dx)^), (3.5) 

of which the first gives rise to a trilinear interaction. As 
these interactions are derivative their largest effects will 
be seen in the high px regime, which we exploited in [l^ 
in order to suppress backgrounds to a manageable level. 
If we consider a strongly interacting SU{N) gauge the- 
ory, then there will be iV^ — 1 gauge fields, and the the- 
ory will be approximately conformal if there are ^ 11 A'^ 
fiavors of Weyl fermion. Taking A^ = 4, 5, 6 we obtain 
auv = 0.033,0.053 and 0.076, using results in We 
will initially take Aa = 0.05, but also consider a 'large' 
anomaly coefficient scenario, where we take Aa = 0.2. 

We summarize the parameter values we use regarding 
double dilaton production in Table HI and show c^xGG in 
brackets to indicate that we usually use the value derived 
by matching with the effective field theory, but sometimes 
switch its effect off altogether. 



Parameter 


Value 


Parameter 


Value 


/ 


850 GeV 


A 


3 (SM) 


cs 


7 (4/3) 


cem 


-17/9 (-1.2) 


Aa 


0.05 (0.2) 


CxxGG 


(0) 



TABLE I: Parameters used in the calculation of double- 
dilaton production in Section Till Al 

Figure [5] shows the differential distribution of fi x BR 
for a number of final states, normalized to those of the 
SM; in both the low and high anomaly coefficient scenar- 
ios. The lower panels show the fully composite SM and 
the upper ones the four family scenario. The effects of 
the higher dimensional operators changing the pt spec- 
trum can be seen entering at around 150 GeV. In the 
fully composite case, while the cross-section for those fi- 
nal states involving either 2 or 4 gluons are boosted with 



respect to the SM, the final states that have proved useful 
in previous dihiggs analyses are suppressed relative to the 
SM, even though the total cross-section for xx is consid- 
erably higher. This is due to the double suppression com- 
ing from the factor / p associated with massive final 
states. Although the 77jj final state cross-section is ten 
times the SM rate, the leading order background is still 
too large to make an effective analysis. As it will never 
be feasible to pick out the relatively few gggg or bbgg 
events from the enormous QCD background, one does 
not expect any signal for this particular scenario. One 
possible exception is in the very boosted regime where 
Pt.x ^ 350 GeV, if the effects of higher dimensional op- 
erators are large. 

On the other hand, the suppression factor into massive 
states is smaller in the four family case, and the over- 
all branching ratios are more similar to their SM values. 
While the extra top-partners enhance the total rate, the 
branching ratio to gluons is not so enhanced so as to ren- 
der an analysis impossible. On the contrary, a x BR for 
hhrr and bbW^W~ is approximately an order of magni- 
tude larger than in the SM, a factor which is enhanced 
even more in the high px tail of the distribution. 

In Fig. [7] we show the effects on the px differential dis- 
tribution of varying the anomaly coefficient Aa and the 
dimension 6 coefficient CxxGGy relative to the 'standard' 
case with Aa = 0.05. The yellow line includes only the 
anomalous derivative couplings which appeared in the 
proof of the a-thcorem. Its effect is similarly boosted 
in the low px region where there is lack of destructive 
interference due to the absence of extra box-diagrams. 
We see that the effects of these interactions becomes im- 
portant for pt,x ~ 350 GeV, where it can change the 
cross-section by a factor of a few. The prospects for us- 
ing the di-dilaton final state to constrain the properties 
of the theory's UV completion are thus promising. 

Summary: The cross-section for di-dilaton production 
is much larger than in the Standard Model. However, the 
future LHC prospects for this scenario exhibit a strong 
dependence on ones assumptions about the UV proper- 
ties of the theory. In the fully composite SM, when the 
suppression associated with non-gluonic final states is 
taken into account, all possibly observable final states are 
too suppressed by their branching ratios to give a signal 
at the LHC. On the other hand, in the four-family sce- 
nario the prospects are excellent, with the cross-section 
for reconstructible final states enhanced by up to an order 
of magnitude. This is large enough that one may begin to 
constrain further facets of the UV theory which manifest 
themselves through higher dimensional operators. 



B. Composite di-Higgs production 

The other possibility to have a light SM-like Higgs bo- 
son that we discuss in this work is the composite Higgs 
scenario. The composite Higgs [13, HO] relies on gang- 
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FIG. 6: Comparison of a{xx) x BR(xi)BR(x2) 
to the values of the SM as a function of Pt,x 
for Aa = 0.05 (left panel) and Aa = 0.2 (right 
panel) and cs = 7,CxxGG = 1- The comparison of 
Aa = 0.05, 0.2 is depicted in Fig. [7] 



ing the electroweak interactions as a subgroup of a larger 
spontaneously broken global symmetry group, e.g. 

50(5) ^ 50(4) ~ SU{2)l x SU{2)r , (3.6) 



10 ^ 




PT,x [GeV] 

FIG. 7: Comparison of cr(xx) for different values of Aa and 
CxxGG as a function of pt,x for eg = 7, / = 850 GeV fixed. 
The blue dotted line gives a comparison of Aa = 0.1 to Aa = 
0.05 for cs, CxxGG fixed, which and shows the dependence on 
the trilinear coupling. 



which contains the gauged SU{2)l. Gauging a subgroup 
is tantamount to explicit breaking of the global sym- 
metry, and the (uneaten) Nambu-Goldstone bosons that 
arise from global symmetry brea king pick up a mass from 
a Coleman- Weinberg potential l6l| that involves both 
gauge and fermion loops and breaks electroweak sym- 
metry [63,[6l[6|. 




FIG. 8: Schematic representation of the 2h + ng irreducible 
one- loop (sub) amplitude and for the involved fermion flavors 
in MCHM5, the gluon lines should be understood as off-shell 
currents contributing to e.g. qq hhg. The amplitudes 
involving the trilinear Higgs vertex (i.e. the irreducible h -\- 
ng (sub) amplitudes) are flavor diagonal due to diagonality 
at the gluon vertices AJifj oc 5ij. We include all partonic 
subprocesses in our calculation. 
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FIG. 9: Branching ratios for the mh = 125 GeV Higgs as a 
function of x in MCHM5. 

To incorporate proper hypercharges we need to extend 
the symmetry group to 5(9(5) x U{l)xi and we identify 
hypercharge asY = X + T^ like in other models of strong 
symmetry breaking [ssf . This mechanism is elegantly de- 
scribed by holographic approaches [1^ , where symmetry 
breaking is realized via the Hosotani mechanism [63 | in 
gauge-Higgs unified models. 

The crucial parameter that measures deviations of the 
physical Higgs' couplings to SM matter and parametrizes 
the model's oblique corrections, is given by ^ = w^//^, 
where / is the analogue to the pion decay constant. Con- 
sistency with experimental data can be achieved without 
tuning, which makes this model class a promising candi- 
date for a BSM Higgs sector. In these composite Higgs 
models one generates fermion masses (at least partially) 
via linear mixings with composite fermionic operators in- 
stead of Technicolor-type interactions to avoid bounds 
^ ^ 1. In total, this amounts to a highly modified di- 
Higgs phenomenology compared to the SM expectation, 
which has alread y be en discussed in Refs. [65l - l68l | in some 
detail. In Ref. [631, the effects of the light additional 
fermionic degrees of freedom in the minimal composite 
Higgs model based on Eq. (|3.6|) (referred to as MCHM5) 
have been included to inclusive di-Higgs predictions be- 
yond LET (see also Ref. [Z^l)- The additional fermions 
that run in the gluon fusion loops strongly enhance the 
cross section, and, therefore, can be highly constrained 
by applying the strategies that involve jet recoils in di- 
Higgs production discussed in our previous paper [l^ as 
we will see below. 

MCHM5 introduces a set of composite vector-like 
fermions that form a complete 5 under S'0(5). The 5 
decomposes under the unbroken SU{2)l x SU[2)r, ip = 
52/3 = (2, 2)2/3 + (1, 1)2/3- Obviously, the 52/3 contains 
a weak doublet of fields with the same quantum numbers 
as the left-handed SM quark doublet ql = {tL,bL)^ and 



right-handed top quark, and we can interpret the large 
mass of the top quark as a mixing effect, 

- Cm^ yf{-ipL^'^){^'lpR) + JTlQl/'Ll/'i? 

-f- AlqlQr + AnfLtR + h.c. , (3.7a) 

where the non-linear Higgs field S is parametrized via the 
SO{5)/SO{4:) coset space generators and can be chosen 
(see e.g. Ref. ^) 

E = (0, 0, sin(V/), 0, cos(V/)) • (3.7b) 

Expanding the non-linear sigma model we recover the 
interactions with electroweak gauge bosons as well as the 
Higgs self-couplings relevant to this study 

+ 2!Z!.,„=(^) („.;„.- + j^z„z'-), (3.8) 

i.e. we need to rescale the SM trilinear hVV vertices 
by a factor y/1 — <^ and we have sin^((/i) //) = v. The 
Higgs branching ratios of MCHM5 are depicted in Fig. |9l 
Following Ref. [6^, we do not include another 5_i/3 
multiplet for generating the bottom quark mass, but in- 
clude it by breaking partial compositeness with an ex- 
plicit coupling of the Yukawa-like interactions. Expand- 
ing Eq. p.7p in the mass diagonal basis, we obtain the 
masses of the fermionic mass spectrum and interactions 
hfifj and hh fifj (where i,j run over the heavy fermion 
flavors) which are relevant for di- Higgs (-|- jet) production 
from gluon fusion, which is the dominant production 
mechanism^ . 

In general, the composite Higgs interactions Eq. (|3.7p 
will not be flavor-diagonal in the space of states that con- 
tains the composite multiplet augmented by tL,R, and 
constraints from both direct detection of flavor measure- 
ments are eminent. For the remainder of this section we 
will choose parameter points that are in agreement with 
these constraints to discuss the composite Higgs model's 
implications on di-Higgs and di-Higgs+jet phenomenol- 
ogy following Ref. |69j . 

We take into account all non-diagonal couplings and 
keep the full mass dependence in the calculation beyond 
any approximation. This results in computationally in- 
tense calculations, especially for the pentagon part in 
gg — > ghh and box gg hh (sub)amplitudes where non- 
diagonality of the hfifj vertices increases the feynman 
graph combinatorics. Fig. [51 



^Di-Higgs production from weak boson fusion Fnl is suppressed, 
also because in addition to the hVV vertices the hhVV vertices 
arc rescaled by 1 — 2f with respect to the SM. Unitarization of 
the Vl Vl — > VlVl, qq amplitudes is partially taken over by the 
exchange of techni-p like resonances. These can be studied in the 
weak boson fusion channels [72ll . 
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pp hh bbW+W- 

FIG. 10: Comparison of composite di-Higgs production pt,h 
spectra with the SM for ^ = 0.25. 

The result in comparison to the SM is shown in Fig. [TU] 
for pp — > hh + X production. For a typical mass 
spectrum mt ~ 174 GeV and the lightest composite 
fermion mnghtcst — 1500 GeV we find agreement with 
the enhanced cross sections as reported in Ref. [69j . 
a{hh)/a^^{hh) ~ 3. The phase space dependence of 
this enhancement is rich and non-trivial as a consequence 
of the non-diagonal couplings and additional mass scales 
that show up in the box contributions, which also in- 
terfere with modified trilinear interactions. Hence, it is 
difficult to comment on quantitative similarities of the 
composite Higgs phenomenology for different parameter 
choices. 

However, on a qualitative level, since the composite 
scale needs typically to be large in order to have agree- 
ment with direct searches and flavor bounds, the inclu- 
sive pp ^ hh -\~ X composite phenomenology will be 
dominated by modifications with respect to the SM at 
medium ptm — 100 GeV. This phase space region is 
mostly sensitive to modifications of the tth coupling and 
the modified trilinear h vertex. At large PT,h we ob- 
serve an enhancement due to the presence of new mas- 
sive fermions in the box contributions of the 9 <7-initiated 
subprocesses, which access the protons' valence quark 
distribution. We note that higher order QCD corrections 
are likely to further enhance the cross section prediction 
beyond the naive SM-rescaling [13, [ill ■ 

We find an even larger enhancement of leading order 
pp ^ hh + jet production cross section, with pxj > 
80 GeV 

cr(/i/i + j) ~ 13.0 fb , (3.9) 

for both scenarios shown in Fig. [TU] This result needs to 
be compared to the corresponding LO prediction in the 
SM which is = 2.8 fb, and amounts to an enhance- 
ment of a factor of 4.6. For the fully hadronized search of 
Ref. [l^ this amounts to S/B ^ 7, which is well beyond 



systematic background uncertainties for high luminosity 
searches. 

The relatively larger increase of the one jet-inclusive 
cross section can be understood along the following lines. 
The additional top partners introduce a new mass scale 
to the one-loop amplitude. At large transverse momen- 
tum, the cross section is dominated by continuum hh 
production which mostly proceeds via box diagrams in 
addition to initial radiation. The latter is increased as a 
result of the newly introduced mass scale in comparison 
to the SM, and initial state radiation allows the initial 
state partons to access the large valence quark parton 
distributions. This effect is also visible in the NLO pre- 
dictions of pp ^ hh + X in composite models employing 
the effective theory approximation [73j . 

Summary: 

The composite Higgs scenario is a well-motivated 
model of electroweak symmetry breaking that is consis- 
tent with current flavor constraints and direct searches 
for heavy top partners. Furthermore, composite Higgs 
models typically predict a large enhanced di-Higgs cross 
section, which is further enhanced in for hh + jet final 
state by the introducing a new mass scale to the phe- 
nomenology. While small di-Higgs(-l-jet) rates in the con- 
text of the SM might hinder a determination of the SM 
Higgs potential in case no further indications of physics 
beyond the SM become available, composite di-Higgs 
production will overcome this shortcoming due to its 
large production cross section. Consequently, also for 
extremely heavy top partners, di-Higgs (-f-jet) production 
is going to provide a powerful test of Higgs compositencss 
at the LHC. 



IV. CONCLUSIONS 

A precise determination of the realization of Higgs 
mechanism sui generis is an important task that has to be 
pursued at the LHC, especially after the recent discovery 
of an SM Higgs-like particle. While measurements based 
on single Higgs boson production provide only indirect 
constraints on the realization of electroweak symmetry 
breaking, the partial experimental reconstruction of the 
Higgs potential is indispensable to gain a fuller under- 
standing at a more fundamental level. 

In this paper, we have investigated di-Higgs and di- 
Higgs+jet production in a variety of model classes, whose 
single Higgs production characteristics can account for 
the observation of the new particle at the LHC. Rather 
than employing an agnostic field theory approach" we 
have picked well-motivated examples of realistic BSM 
(scalar) sectors, supplemented by the required fermionic 



II See Ref. [74ll for related discussions. 
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particle content, which generalize the SM Higgs sector in 
two fundamentally distinct ways. 

The first option deals with models with extended Higgs 
sectors predicting new resonant structures in di-Higgs 
production due to the model's two-Higgs doublet char- 
acter. Furthermore, new kinematical configurations can 
provide extra analytical handles in the production of a 
heavy Higgs boson in addition to the light Higgs. In 
portal-inspired scenarios, the determination of the in- 
volved trilinear couplings is important to reconstruct the 
full extended portal potential for parameter points where 
the two Higgs bosons are not too widely separated in 
mass. In the MSSM, a corresponding measurement fa- 
cilitates the reconstruction of the Higgs sector mixing 
angles a and /3, and hence provides indirect constraints 
on the stop masses and mixing parameter At- This can 
be achieved by separating the resonant contribution from 
continuum production via invariant mass cuts, and ap- 
plying boosted [l^ and unboosted [l3| analysis strategies 
to the different samples. 

The resonant models are contrasted to realizations of 
the Higgs mechanism where the "Higgs" boson arises as 
a pseudo Nambu-Goldstone mode of some spontaneously 
broken symmetry. The agreement of current observa- 
tions with the SM Higgs predictions requires the pseudo- 
Nambu-Goldstone boson to have similar couplings as the 
SM Higgs boson. Along with composite Higgs models 
this leaves only the pscudo-dilaton as a second option. 

The former case implies interpreting the entire Higgs 
doublet as a set of Nambu-Goldstone fields. Realistic 
composite Higgs scenarios predict strongly-enhanced di- 
Higgs and di-Higgs-l-jet cross sections. 

In models with an approximate conformal invariance, 
symmetry breaking can be triggered at scales consider- 



ably higher and spontaneous breaking of conformal in- 
variance introduces a new light state to the low energy 
effective theory, which has similar properties as the SM 
Higgs boson as a consequence of Eq. p.ip : the pseudo- 
dilaton. For both the composite and the dilaton option, 
there are parameter choices such that current observa- 
tions can be accounted for. It is their highly modified 
di-Higgs phenomenology which can effectively discrimi- 
nates these possibilities depending on the further partic- 
ularities of the conformal sector, and facilitates an LHC 
measurement of the involved couplings and parameters 
in case of the composite Higgs model. Pseudo-di-dilaton 
production can be buried in a large hadronic background 
with no kinematic handles to reconstruct the preferred 
dilaton decay to gluons. In this sense, the absence of 
a "traditional" di-Higgs phenomenology could be inter- 
preted as evidence for a dilatonic realization. Interpret- 
ing the presence of a large di-Higgs(-l-jet) production 
cross section is more involved, and could be evidence for 
a fourth-family (or more complicated) realization of the 
pseudo-dilaton model, but may also be consistent with a 
composite Higgs. 

It is clear from our analysis that, no matter what 
governs the dynamics of the newly-discovered boson, its 
multi-production phenomenology, which can be studied 
at the LHC in sufficient detail, will provide a clear im- 
age of its role in the mechanism of electroweak symmetry 
breaking. These findings will further consolidate with an 
LHC luminosity upgrade [20| . 
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